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Abstract 


Two  problems  related  to  the  Royal  Thai  Navy's  recruitment  effort  are 
modeled  and  solved  by  mathematical  programming.  The  first  problem  is  the 
optimal  assignment  of  draftees  to  branch  naval  bases,  which  is  formulated 
as  a  transportation  problem.  The  second  problem  is  the  optimal  routing 
of  ships  for  transporting  draftees  from  the  branch  bases  to  the  main 
base,  which  is  solved  with  an  exact  integer  programming  formulation. 


S- 


Introduction 


There  is  consuls ory  military  service  for  men  in  Thailand.  Inductions 
take  place  four  times  a  year,  with  the  men  who  live  in  the  coastal  provinces 
assigned  to  the  Navy. 

When  a  man  is  called  into  the  Thai  Navy,  he  first  reports  to  a  drafting 
center  in  his  home  locality.  He  is  then  transported  by  land  vehicle  to  a 
naval  base.  If  he  is  from  a  northern  province,  this  base  is  the  main  naval 
base  in  Sataheep,  near  Bangkok.  If  he  is  from  a  southern  province,  the  draftee 
is  first  brought  overland  to  a  branch  naval  base  and  then  transported  by  ship 
to  the  main  base.  (See  the  map  of  Figure  1.)  In  the  southern  provinces  there 
are  36  drafting  centers  and  4  branch  bases.  This  paper  addresses  two  problems 
relating  to  the  assignment  and  transportation  of  draftees  from  the  southern 
provinces : 

Problem  1:  How  many  men  from  each  drafting  center  should  be  assigned  and 
transported  to  each  branch  base? 

Problem  2:  Given  an  available  fleet  of  ships  stationed  at  the  main  base, 
which  of  those  ships  should  be  used  and  how  should  they  be  routed  so  as 
to  pick  up  all  the  men  at  the  branch  bases  and  transport  them  to  the  main 
base? 

The  second  problem  is  the  more  difficult. 

We  shall  present  a  linear  progranming  model  which  is  used  for  solving  the 
first  problem  and  an  integer  programming  model  which  is  used  for  solving  the 
second.  The  two  problems  are  not  independent  because  the  solution  to  Problem 
1  determines  the  transport  requirements  for  Problem  2. 


mould  be  possible  to  build  a  unified  model  to  solve  both  problems  simultane¬ 
ously,  this  would  be  inappropriate  under  the  Thai  Navy's  current  organiza¬ 
tional  structure.  The  personnel  decisions  of  Problem  1  and  the  ship  routing 
decisions  of  Problem  2  are  made  by  different  people  in  different  places  at 
different  times. 

Problem  1:  Assigning  Draftees  to  Branch  Bases 

Problem  1  is  modeled  as  a  standard  transportation  problem. 

We  are  given: 

bj  =  the  number  of  men  to  be  transported  from  drafting  center  j, 
aj  =  the  capacity  (number  of  men)  of  branch  base  i, 

Cij  =  cost  per  man  for  transport  to  base  i  from  center  j, 
and  we  must  determine 

xij  =  the  number  of  men  to  be  transported  to  base  i  from  center  j. 

The  model  is 

°ijXij 

Si  t« 

Exij  <  a*,  all  i 

£xij  =  bj>  a11  i 
i 

with  all  Xij  >  o.  The  vehicles  used  for  transporting  the  men  are  sufficiently 

small  so  that  it  is  not  necessary  to  use  a  more  elaborate  model  (e.g,  a  fixed- 

charge  transportation  model)  for  Problem  1.  The  results  of  Problem  1  which 

are  inputs  to  Problem  2  are 

<*i  =  Z*ij 

j 

which  is  the  number  of  men  assigned  to  branch  base  i.  The  data  and  numerical 
results  for  a  specific  instance  of  Problem  1  (Summer,  1983)  are  given  in 


Tables  1  and  2 
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TABLE  1.  DATA  FOR  TRANSPORTATION . PROBLEM 


Drafting 

Center 

Distances 

to  Branch  Bases 

Number  of  Men 
at  Center 

B1 

B2 

B3 

B4 

PC 

25 

100 

170 

270 

98 

KR 

30 

85 

165 

260 

82 

PN 

10 

87 

150 

250 

123 

SW 

23 

77 

140 

240 

54 

RN 

50 

75 

150 

240 

118 

CY 

77 

20 

90 

180 

98 

TK 

112 

29 

87 

175 

73 

TP 

130 

70 

120 

200 

66 

BD 

100 

0 

75 

170 

102 

SC 

125 

50 

40 

140 

57 

TS 

130 

60 

23 

125 

61 

PG 

160 

77 

112 

175 

104 

KL 

170 

90 

120 

178 

60 

TL 

180 

110 

120 

170 

58 

PK 

195 

120 

125 

70 

98 

KB 

175 

82 

80 

138 

107 

CW 

150 

50 

28 

123 

92 

TG 

170 

75 

25 

100 

107 

RB 

175 

80 

20 

87 

93 

KM 

180 

80 

48 

85 

62 

RT 

210 

125 

73 

50 

69 

PL 

215 

125 

70 

50 

88 

TR 

210 

120 

75 

77 

93 

KT 

220 

127 

82 

78 

86 

SN 

227 

135 

74 

28 

74 

HY 

267 

175 

110 

12 

103 

PR 

260 

167 

112 

70 

52 

ST 

282 

185 

132 

55 

77 

TA 

282 

208 

133 

37 

59 

SD 

280 

190 

130 

31 

68 

PT 

302 

220 

150 

60 

83 

YR 

310 

228 

162 

73 

61 

SB 

325 

235 

175 

81 

82 

NW 

348 

270 

200 

107 

% 

YL 

325 

230 

168 

70 

72 

BT 

360 

275 

203 

100 

53 

Base 

Capacities  ICOO 


1100 


1300 


1500 


TABLE  2.  OPTIMAL  SOLUTION  FOR  TRANSPORTATION  PROBLEM 


To 

Number  of  Men 

Cost/Man 

Extension 

B1 

98 

25 

2450 

B1 

82 

30 

2460 

B1 

123 

10 

1230 

B1 

54 

23 

1242 

B1 

118 

50 

5900 

B2 

98 

20 

1960 

B2 

73 

29  . 

2117 

B2 

66 

70 

4620 

B2 

102 

0 

0 

B2 

104 

77 

8008 

B2 

60 

90 

5400 

B2 

58 

110 

6380 

B2 

98 

120 

11760 

B3 

57 

40 

2280 

B3 

61 

23 

1403 

B3 

107 

80 

8560 

B3 

92 

28 

2576 

B3 

107 

25 

2675 

B3 

93 

20 

1860 

B3 

62 

48 

2976 

B3 

93 

75 

6975 

B4 

69 

50 

3450 

B4 

88 

50 

4400 

B4 

86 

78 

6708 

B4 

74 

28 

2072 

B4 

103 

12 

1236 

B4 

52 

70 

3640 

B4 

77 

55 

4235 

B4 

59 

37 

2183 

B4 

68 

31 

2108 

B4 

83 

60 

4980 

B4 

61 

73 

4453 

B4 

82 

81 

6642 

B4 

96 

107 

10272 

B4 

72 

70 

5040 

B4 

53 

100 

5300 

Total  Cost  =  149,551 


Number  of  Men  Assigned 


Problem  2:  Routing  Ships  for  Transporting  Draftees  to  the  Main  Base 

The  problem  of  optimally  deploying  the  available  fleet  of  ships  to  transport 
draftees  to  the  main  base  is  a  type  of  vehicle  routing  problem.  There  is  a  large 
number  of  available  approaches  to  these  problems.  (See,  e.g.,  Bell  et.  al. 
(19831  and  Bodin  et.  al.  (19831).  The  model  given  below  is  related  to  the  set- 
partitioning  approach,  in  that  we  first  create  a  list  of  possible  routes  that 
ships  can  take  and  then  model  the  problem  in  terms  of  variables  representing  the 
number  of  times  each  route  is  used.  Since  ships  are  not  identical,  these  var¬ 
iables  must  also  denote  which  type  of  ship  is  chosen. 

In  most  vehicle  routing  problems,  the  objective  function  is  the  total  dis¬ 
tance  traveled  or  some  concomitant  cost  measure.  The  Thai  Navy's  situatior 
different.  Their  primary  objective  is  to  minimize  the  number  of  ships  req  id 
for  transporting  all  the  draftees.  This  is  because  defending  the  country,  ’* 
transporting  draftees,  is  the  principal  mission  of  the  Navy  fleet.  As  few  ves¬ 
sels  as  possible  should  be  assigned  to  this  duty,  so  that  the  remaining  ships  can 
be  available  for  defense.  Since  we  anticipate  that  numerous  solutions  would  al¬ 
ternatively  optimize  this  primary  objective,  we  use  distance  minimization  as  a 
secondary,  tie-breaking  objective. 

The  fleet  available  for  draftee  transport  at  the  time  of  this  study  con¬ 
sisted  of  9  vessels  of  3  different  classes.  This  fleet  is  described  in 
Table  3. 

TABLE  3.  FLEET  DESCRIPTION 

Class  Capacity  of  Ship  Number  of  Bases  that  Class 

of  Ship  (Number  of  Men)  Ships  in  Class  Can  Visit 


From  this  information  we  generated  23  possible  voyages  for  the  model  to 
consider  using.  A  voyage  is  defined  as  a  particular  class  of  ship  assigned  to 
sail  a  particular  route.  The  23  voyages  are  described  in  Table  4.  Node  0  in 


the  route  specifications  means  the  main  base. 


TABLE  4.  VOYAGE  SPECIFICATION 


Voyage 

Class 

Route 

Distance  (Ion) 

Capacity  (men) 

1 

1 

0-1-0 

370 

100 

2 

2 

0-1-0 

370 

200 

3 

2 

0-2-0 

460 

200 

4 

2 

0-3-0 

600 

200 

5 

2 

0-1-2-0 

515 

200 

6 

2 

0-1-3-0 

640 

200 

7 

2 

0-2 -3-0 

665 

200 

8 

2 

0-1-2-3-0 

720 

200 

9 

3 

0-1-0 

370 

600 

10 

3 

0-2-0 

460 

600 

11 

3 

0-3-0 

600 

600 

12 

3 

0-4-0 

750 

600 

13 

3 

0-1-2-0 

515 

600 

14 

3 

0-1-3-0 

640 

600 

15 

3 

0-1-4-0 

810 

600 

16 

3 

0-2-3-0 

665 

600 

17 

3 

0-2-4-0 

805 

600 

18 

3 

0-3-4-0 

800 

600 

19 

3 

0-1-2-3-0 

720 

600 

20 

3 

0- 1-2-4-0 

860 

600 

21 

3 

0-1-3-4-0 

840 

600 

22 

3 

0-2 -3-4-0 

865 

600 

23 

3 

0-1-2-3-4-0 

920 

600 

The 

method  used  for  generating  the 

set 

of  available  voyages  is  the  fol- 

lowing. 

For  each  ship 

class,  we  enumerate  every  subset  of 

bases  that  can  be 

visited  i 

an  a  voyage  and 

then  we  compute 

the 

shortest  route 

for  that  subset . 

All  voyages  worth  considering  are  thereby  generated.  If  the  ship  routing  prob- 

lem  were 

posed  by  another  fleet  whose  number 

of  classes  and 

number  of  bases 

were  much  greater  than  the  Thai's,  then  this  exact  approach  might  become  com¬ 
putationally  infeasible. 


The  variables  in  our  integer  programming  model  for  Problem  2  are: 


zj  =  number  of  times  voyage  j  is  used,  and 

yij  =  number  of  men  transported  from  base  i  to  the  main  base  via 
voyage  j. 

The  data  for  the  integer  programming  model  are: 

=  set  of  voyages  that  visit  base  i, 

Sfc  =  set  of  voyages  that  use  a  ship  of  class  k, 

Bj  =  set  of  bases  visited  on  voyage  j, 
dj  =  number  of  men  to  be  transported  from  base  i, 
nfc  ■  number  of  ships  available  of  class  k, 
cj  =  capacity  of  voyage  j,  and 

rj  =  distance  traveled  by  voyage  j  in  tens  of  thousands  of  kilometers. 
The  model  is 


min  Z  (1  ♦  r. 

+  J  J 

s.t. 


jivi  =  “i-  a11  1 

tfej  7ij  ‘  °jV  511  J 

jisk  zj  ‘  a11  k 


Z 


j 


>  0  and  integer 


>  0. 


The  objective  function,  as  noted,  primarily  serves  to  minimize  the  number 
of  voyages  required.  The  rj  have  small  values  and  serve  only  to  break  ties. 
The  first  constraint  set  insures  that  all  the  drafted  men  at  base  i  are  trans¬ 
ported,  using  voyages  that  in  fact  pass  trough  base  i.  The  second  constraint 
set  insures  that  the  number  of  men  transported  via  voyages  of 


type  j  is  within  the  capacity  that  will  be  made  available.  The  third  con¬ 
straint  set  insures  that  each  ship  In  the  fleet  is  used  at  most  cnce.  A  ship 
can  not  be  reused  because  all  the  draftees  have  to  start  their  training  at  the 
main  base  at  the  same  time. 

Software  Employed 

The  transportation  model  of  Problem  1  is  of  course  a  standard  problem  for 
which  numerous  effective  packages  exist.  In  our  case  we  used  GNET,  the  ca¬ 
pacitated  transshipment  problem  solver  of  Bradley,  Brown  and  Graves  [19771 .  We 
wrote  a  problem  generator  and  report  writer  in  FORTRAN  to  use  in  conjunction 
with  this  program,  and  we  used  the  DEC-10  computer's  T0PS-10  operating  system 
for  managing  the  necessary  file  interfaces. 

The  integer  programming  model  in  the  particular  instance  we  solved  was 
small  enough  to  handle  with  the  general  purpose  solver  UNDO  of  Schrage 
[1981].  The  model  formulation  above  for  the  Thai  fleet  of  Tables  3  and  4  has 
23  integer  variables  zj  and  45  continuous  variables  y^j.  The  zj  variables  can 
take  cn  any  integer  value  up  to  and  including  the  number  of  ships  in  the  class 
assigned  to  voyage  j.  UNDO,  like  many  other  commercial  integer  programming 
solvers,  requires  that  all  the  integer  variables  be  binary.  Therefore,  we  had 
to  reformulate  the  model  with  the  zj  replaced  by  the  binary  expansions: 

Zj  3  ZJ1  *  2zj2,  J  s  8, 

Zj  3  zjl  +  2zJ2  ♦  4zj3,  j  *  9,..., 23 

This  resulted  in  a  model  with  61  binary  variables. 

Another  problem  generator  was  written  to  automatically  generate  the 


voyages  and  create  this  input,  so  that  other  instances  of  the  model  can  be 
treated  routinely.  The  results  for  the  model  instance  based  on  Tabls  3  and 
4  are  given  in  Tables  5  and  6. 


TABLE  5.  OPTIMAL  SOLUTION  TO  SHIP  ROUTING  MODEL 


Nonzero  Variable 

Value 

Interpretation 

Z{2,2) 

1 

Use  voyage  2  twice. 

2(5,1) 

1 

Use  voyage  5  once. 

Zd0,l) 

1 

Use  voyage  10  once. 

2(11,1) 

1 

Use  voyage  11  once. 

Z ( 12 , 1 ) 

1 

Use  voyage  12  once. 

Z(18,l) 

1 

Use  voyage  18  once. 

Y(l, 2) 

400 

Transport  200  men  cn  each 

voyage  2. 

Y(l, 5) 

75 

Transport  75  men  from  B1  on  voyage  5. 

Y(2,5) 

59 

Transport  59  men  from  B2  on  voyage  5. 

Y(2,10) 

600 

Transport  600  men  from  B2 

on  voyage  10. 

Y(3,ll) 

595 

Transport  595  men  from  B3 

on  voyage  11. 

Y(3,18) 

77 

Transport  77  men  from  B3  on  voyage  18. 

Y(4,12) 

600 

Transport  600  men  from  B4 

on  voyage  12. 

Y ( 4,18) 

523 

Transport  523  men  from  B4 

on  voyage  18. 

TABLE  6. 

OPTIMAL 

SHIPPING 

SCHEDULE 

Class 

Voyage 

of  Ship 

Route 

Men 

Carried 

Distance 

2 

2 

0-1-0 

200 

370 

2 

2 

0-1-0 

200 

370 

5 

2 

0-1-2-0 

75+59 

515 

10 

3 

0-2-0 

600 

460 

11 

3 

0-3-0 

595 

600 

12 

3 

0-4-0 

600 

750 

18 

3 

0- 3-4-0 

77+523 

800 

Total  number 

of  voyages 

=  7 

Total  distance 

=  3865  km 

Discussion  of  Results 

According  to  Thai  Naval  authorities,  the  results  obtained  by  the  math¬ 
ematical  programming  models  are  superior  to  the  results  obtained  through  the 
manual  procedures  currently  employed.  In  the  instance  of  the  transpor-tatian 
model  reported  above,  it  turned  out  that  just  sending  each  man  to  the  nearest 
base,  as  would  be  done  by  the  manual  procedure,  turned  out  to  be 
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feasible.  However,  in  other  instances,  this  procedure  has  resulted  in  vio¬ 
lations  of  the  base  capacities,  which  necessitated  additional  transportion  of 
men  from  oversubscribed  bases  to  undersubscribed  bases.  The  authorities  were 


pleased  to  see  that  the  cost  and  time  delays  of  these  extra  trips  could  be 
avoided. 

With  respect  to  the  ship  routing  model,  again  the  authorities  were  pleas¬ 
ed  to  see  our  recommendation  of  an  efficient  deployment  of  their  fleet.  Data 
gathering  is  currently  in  process  to  determine  how  many  voyages  and  kilometers 
could  have  been  saved  if  the  mathematical  programming  models  had  been  avail¬ 
able  for  the  past  few  years. 

Conclusion 


Two  mathematical  programming  models,  a  standard  transportation  model  and 
a  new  and  exact  integer  programming  formulation,  have  been  shown  to  be  ef¬ 
fective  for  optimally  routing  ships  and  assigning  personnel  for  Thai  Naval 
recruitment.  The  models  were  solved  by  general  purpose  software.  Our  ap¬ 
proach  is  potentially  applicable  to  other  situations,  however  it  is  likely 
that  special  purpose  software  would  be  needed  for  the  integer  programming  pro¬ 
blem  if  larger  fleets  are  considered. 
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